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We propose a scheme for generating two-dimensional turbulence in harmonically trapped atomic 
condensates with the novelty of controlling the polarization (net rotation) of the turbulence. Our 
scheme is based on an initial giant (multicharged) vortex which induces a large-scale circular flow. 
Two thin obstacles, created by blue-detuned laser beams, speed up the decay of the giant vortex 
into many singly-quantized vortices of the same circulation; at the same time, vortex-antivortex 
pairs are created by the decaying circular flow past the obstacles. Rotation of the obstacles against 
the circular flow controls the relative proportion of positive and negative vortices, from the limit 
of strongly anisotropic turbulence (almost all vortices having the same sign) to that of isotropic 
turbulence (equal number of vortices and antivortices). Using the new scheme, we numerically 
study quantum turbulence and report on its decay as a function of the polarization. We finally 
present a phenomenological model for the decay rate of vortex number which hts our numerical 
experiment curves, with the novelty of taking into account polarization time-dependence. 

PACS numbers: 03.75.Lm, 67.25.dk, 67.85.De 


I. INTRODUCTION 

The study of quantum turbulence is heavily motivated 
by liquid helium (^He and ^He) experiments [HE]- A 
striking discovery has been that, under appropriate forc¬ 
ing, quasi-classical behavior arises displaying statistical 
properties characteristic of ordinary turbulence; an ex¬ 
ample is the celebrated Kolmogorov —5/3 scaling of the 
energy spectrum [3] which suggests the existence of a 
classical energy cascade from large to small length scales. 
Under other conditions, a different kind of turbulence 
(called ‘ultra-quantum turbulence’ or ‘Vinen turbulence’) 
has also been found m, characterized by random tan¬ 
gles of vortices without large-scale, energy-containing 
flow structures. Quantum turbulence experiments are 
also performed in atomic Bose-Einstein condensates [i- 
[8] ; the relative small size of these condensates (compared 
to flows of liquid helium or of ordinary fluids) limits the 
study of scaling laws but offers opportunities to study 
minimal processes that also take place in larger systems 
(e.g. vortex interactions, vortex reconnections, vortex 
clustering) with greater experimental controllability and 
more direct visualization than in liquid helium. 

Atomic condensates are also ideal systems to study 
two-dimensional (2D) turbulence [9], a problem with im¬ 
portant applications to oceans, planetary atmospheres 
and astrophysics. In these systems, reduced dimension¬ 
ality may arise from strong anisotropy, stratification or 
rotation (via the Taylor-Proudman theorem). From the 
physicist’s point of view, the dynamics of 2D turbulence 
is very different from 3D m- The existence (besides 
the kinetic energy) of a second inviscid quadratic invari¬ 
ant - the enstrophy - implies that a downscale enstrophy 
transfer is accompanied by an upscale energy cascade; in 


other words, in 2D turbulent flows the energy flows from 
small to large length scales rather than vice-versa as in 
3D turbulence. With the possible exception of soap films 
m, 2D flows which can be created in the laboratory are 
only approximations. However, using suitable trapping 
potentials, atomic condensates can be easily shaped so 
that vortex dynamics is 2D rather than 3D. Unlike liq¬ 
uid helium, in atomic condensates 2D quantum vortices 
can be directly imaged, and, unlike classical systems, the 
motion of such 2D vortices is not hindered by viscous 
effects or friction with the substrate. 

Several works have explored the generation of turbu¬ 
lence in 2D condensates. The 2D energy spectrum and 
scaling laws have been computed in numerical simula¬ 
tions [EHIll, and the problem of what should be the 
quantum analogue of the classical enstrophy has been 
raised. In Ref. [8] vortices were nucleated by small-scale 
stirring of a laser spoon, after which a persistent cur¬ 
rent was verified both experimentally and through nu¬ 
merical simulations, suggesting transfer of incompressible 
kinetic energy from small to large length scales. Emer¬ 
gence of large-scale order from vortex turbulence was also 
observed m as predicted by the ‘vortex gas’ theory of 
Onsager. A similar set-up was used to explore vortex 
shedding and annihilation processes in both experiments 
min] and simulations [18]. The effect of stirring laser 
beams with different shapes or along different paths was 
investigated in Refs. [T9H22]- However, in all of cited 
cases, vortices have always been generated in such a way 
that the number of positive and negative vortices is ap¬ 
proximately the same; in other words, all vortex configu¬ 
rations which have been investigated had approximately 
zero polarization. Since irrotational flow is a hallmark 
property of superfluidity, the polarization of the vortex 
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configuration (i.e. the relative proportion of positive and 
negative vortices) plays the role of net rotational angular 
velocity of a classical fluid, so it is important to explore 
its effects on the properties of turbulence. 

In this work, we propose a new scheme for generating 
2D quantum turbulence in atomic condensates. The nov¬ 
elty of our scheme, which is based on a giant vortex as 
the initial state, is control over polarization of the turbu¬ 
lence, which can be interpreted as the classical rotation 
of the entire flow. One of the most important properties 
of turbulence is its decay, because the growth of the tur¬ 
bulence or its character in a steady state may depend on 
how it is forced, whereas the decay is an intrinsic prop¬ 
erty of the dynamics. We shall report the decay of 2D 
quantum turbulence as a function of the polarization. 


II. GIANT VORTEX AND SMALL PINS 

Multicharged vortices with circulations as large as 60 
quanta have already been produced in condensates us¬ 
ing dynamical methods, as consequences of rapid rota¬ 
tions of the confining trap [23] . Another route to achieve 
these highly excited states is using phase-engineering 
techniques, such as those described in Refs [241127] . In 
these cases, quanta of angular momentum are added to 
the condensate by adiabatically inverting the direction of 
the magnetic bias-field which composes the usual loffe- 
Pritchard magnetic trap. Up to this date, only charges 
below 10 quanta were produced using their proposed set¬ 
ups. However, an improvement on the method, known as 
the ‘vortex-pump’, has been described in Refs. [28II3Q] . 
In practical terms, a hexapole magnetic field is super¬ 
posed to the loffe-Pritchard magnetic trap, allowing vor- 
ticity to be cyclically pumped into the condensate, and 
generating giant vortices. Progress in this direction has 
been done in recent experiments with synthetic magnetic 
monopoles m- 

A giant vortex at the center of a harmonically trapped 
condensate can be described by a single-particle wave- 
function of the form '0(r) = where /(r) is the 

wave-function’s amplitude, r = (r, 0, z) is the position 
in cylindrical coordinates, and a large winding number 
K corresponds to a large angular momentum. Such gi¬ 
ant vortices are dynamically unstable [29l |32] , and split 
into singly-quantized vortices. Being parallel to one an¬ 
other, these singly-quantized vortices impose a strongly 
azimuthal flow to the condensate. During the following 
evolution, some vortices of the opposite polarity may be 
generated by occasional large-amplitude density waves, 
but these events are rare, and do not change the main 
property of the flow resulting from the decay of a giant 
vortex configuration: the strong polarization of the vor- 
ticity - almost all vortices have the same sign. 

The scheme that we propose uses blue-detuned lasers 
m to perturb this initial state with two diametrically 
opposite laser beams, creating thin obstacles (which we 
refer as pins) with width a of the order of magnitude of 


the healing length ^ (two pins are enough to homogenize 
the vortex distribution). The pins perturb the initial gi¬ 
ant vortex, accelerating its decay; they also deflect the 
large azimuthal flow, generating vortex-antivortex pairs 
[I81[33H35]. To control the effect of the pins, we move 
them at constant angular velocity uj in the direction op¬ 
posite to the main azimuthal flow. 

III. MODEL 

The dynamics of our system is dictated by the 2D 
Gross-Pitaevskii equation (GPE). We introduce dimen¬ 
sionless variables based on the trapping potential of fre¬ 
quency ujQ, measuring times, distances, and energies in 
units of y/h/mujQ and huuo respectively, where m 

is the mass of one atom and h is the reduced Planck’s 
constant. The resulting dimensionless GPE is 

= + + ( 1 ) 

where the time-dependent wavefunction '0(r, t) is nor¬ 
malized so that J \'ip\‘^d^r = 1. The external poten¬ 
tial is U(r,t) = Utrap(r) + Upins(r,t), where Utrap(r) = 
(x^+^^)/2 and Upins(i*,t) = U+(r, t) + U_(r, t) represent 
respectively the trapping potential which confines the 
condensates and the pins which perturb the initial giant 
vortex. The terms V± (r, t) = Vq exp { — | r — r± (t) p /2cr^ } 
with r±(t) = [±xo cos (cjt), ^0 sin (cjt)] are diametrically- 
opposite, thin, Gaussian potentials of width a = ^ which 
rotate clockwise (against the flow of the initially imposed 
giant vortex) at constant angular velocity uj. The quan¬ 
tity C = 2\/^N{a/ az) parametrizes the two-body col¬ 
lisions between the atoms, where N is the total number 
of atoms, a the scattering length, and az the axial har¬ 
monic oscillator’s length; we choose C = 17300. The 
chemical potential /i is introduced to guarantee normal¬ 
ization of the wave-function, and the amplitude of the 
pins is Vb ~ 1.43/i. In homogeneous systems (U = 0) 
the healing length is found by balancing kinetic and in¬ 
teraction energies terms in the GPE. In a harmonically 
trapped condensate, the healing length can be defined 
with reference to the density at the center of the trapped 
condensate in the absence of any vortex or hole. In our 
dimensionless units, we obtain ^ ^ 0.13, and ttf ^ 74^ 
for the Thomas-Eermi radius. 

Our choice of dimensionless parameters corresponds 
to typical HSl HU experiments with ^^Na condensates 
(scattering length a = 2.75 nm, atom mass m = 3.82 x 
10“^^kg) with N = 1.3 X 10^ atoms, radial and axial 
trapping frequencies c^o = 27r x 9 Hz and ujz = 27r x 400 
Hz, radial and axial harmonic oscillator’s lengths uq = 
y/h/mujQ ~ 7.1 /am and az = ^/h/mujz ~ 1.0 /im, for 
which the dimensional healing length is ^ = 0.13ao ~ 0.9 
/im; the laser beam would have then a Gaussian 1/e^ ra¬ 
dius of iCo = 2(7 ~ 1.8 /im. Blue-detuned Gaussian laser 
beams have been used as pins in a series of experiments 
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with highly-oblate BECs [ITl[22j|36] . Particularly in [36] . 
a laser beam of width wq ^ 2 /im was used to stir a 2D 
^^Rb condensate, similarly to what we propose, main¬ 
taining a circular motion with the help of piezo-driven 
mirrors. 

In order to define our initial state, a circulation of 37 
quanta (i.e. winding number k. = 37) is initially im¬ 
printed around the center of the Thomas-Fermi profile, 
thus imposing an initial counter-clockwise circular flow. 
Changing t into —it in Eq. Q, we shortly evolve the state 
for t = 0.09 in imaginary-time description, guaranteeing 
a fixed phase of 27rK in the center of the condensate and 
adjusting the density to the presence of the pins. We 
then compute the evolution in real time. 

By substituting t in Eq.[^ we are left with a 

phenomenological dissipative GPE (dGPE), where 7 is a 
dissipation constant which models the interaction of the 
condensate with the surrounding thermal cloud. This 
equation can be used to investigate the effect of finite 
temperature in our system. With this aim, we also repeat 
our simulations using dGPE with 7 = 3.0 x 10“^, a typ¬ 
ical value of dissipative parameter [37II39] . particularly 
chosen for the for experimental realistic case found in 
current experiments [ISllIH]. Summarizing beforehand, 
we find the same overall behavior for both dissipation-less 
and this specific dissipative case. 

All numeric simulations are performed in the 2 D do¬ 
main —25 < < 25 on a 512 x 512 grid using the 4th 

order Runge-Kutta method in Fourier space with the help 
of XMDS go]. 


IV. RESULTS 
A. Creating Polarized Flow 

We simulate the real-time evolution of the system for 
different values of the pins’ angular velocity: cj = 0, 7r/16, 
7 r/ 8 , 7 r/ 6 , and 7 r/ 4 . A series of snapshots for the case of 
cj = 0 is shown in FIG. to exemplify a typical run. 
The initial large hole at the center of the figure is the 
core of the giant vortex. The two small holes (north and 
south of the giant hole) are the two stationary pins. The 
critical velocity Vc for the creation of a vortex-antivortex 
pair depends on the barrier’s shape HHun] and also on 
inhomogeneities of the system [T7]. Typically Vc/c ~ 
0.1 — 0.4 for infinitely high cylindrical barriers, where 
c is the local speed of sound. Since our barriers (the 
pins) are either stationary or rotate against the main 
flow, vortex shedding is a dissipative mechanism which 
slows down the superfluid’s azimuthal flow and removes 
angular momentum. 

Besides generating vortices of opposite sign, the pins 
act as perturbation to the giant vortex and quicken its 
decay process; for example, a wave front which perturbs 
the core of the giant vortex is visible at time t = 0.9 
in FIG. The decay of the giant vortex takes place 
via deformation of the core, which becomes elliptical be¬ 



FIG. 1. (Color online). Density plots of the condensate at 
different times t for oj = 0 (non-rotating obstacles). Regions 
of large/low density are displayed in white/black respectively. 
Red triangles and blue circles identify positive-charged and 
negative-charged vortices respectively. The giant vortex de¬ 
cays by injecting a large number of singly-quantized (positive) 
vortices into the condensate, whilst the pins generate vortex 
pairs, as can be clearly seen at time t = 0.9. 


fore vanishing, and injecting a large number of positive, 
singly-quantized vortices into the condensate. At the 
same time, vortex-anti vortex pairs are created by the flow 
past the pins. This process continues until the large az¬ 
imuthal flow is lower than the critical velocity Vc] at that 
point the giant vortex has disappeared, and the pins are 
practically unable to generate further vortices. There¬ 
fore, after this slowdown and due to their small sizes, 
the pins are practically irrelevant to the vortex dynamics 
(apart from occasional creation of pairs in the fast rotat¬ 
ing case, cj = 7r/4). In spite of that, in order to study the 
vortex number decay, we simply remove them at t = 82 
and allow for longer simulations. 

We perform a phase-unwrapping procedure and, by de¬ 
tecting windings of ±27r around small closed paths (pla- 
quettes) on the phase-profile [19] , we are allowed to count 
the numbers and N~ of positive and negative singly- 
quantized vortices in the system (anticlockwise and clock¬ 
wise circulation respectively). This vortex detection al¬ 
gorithm uses a density-cut criterion (^ 0.75 of 7 /^’s mean- 
value) to avoid detection of ghost-vortices. Given the 
initial giant vortex (which is multicharged and therefore 
not detected by our vortex-detecting algorithm), depend¬ 
ing on the value of cj, there can be an imbalance of 
and N~ throughout the evolution. Vortices can be ex¬ 
pelled from the condensate due to their mutual interac- 
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FIG. 2. (Color online), (a) total number of vortices, iVtot 
vs time t; (b) polarization P vs time t. The vertical dashed 
line marks the time (t = 82) we use to make initial states 
for longer simulation without the pins. The curves are dis¬ 
tributed in increasing value of uj from bottom to top, for 
the top plot, and conversely, for the bottom plot. In (a), 
at t 8, from bottom to top, the curves refer respectively 
to Ct; = 0,7r/16,7r/8, tt/G, and 7r/4. In (b), from top to bot¬ 
tom, the curves refer respectively to the same latter series of 
increasing values of cj. 


tion, spiral out of the condensate because of dissipation, 
or undergo vortex-pair annihilation processes. In our par¬ 
ticular finite-temperature simulation, we verify that the 
chosen experimental realistic value of the dissipation pa¬ 
rameter 7 is small enough that, on the time scale ana¬ 
lyzed (and compared to the dissipation-less simulations), 
dissipation-induced spiraling out of individual vortices is 
less important than vortex interactions or annihilations. 

After the decay of the initial giant vortex, the imbal¬ 
ance of positive and negative vortices is measured by the 
polarization P = {N~^ — A^“)/(A^+ + which takes 

maximum/minimum values (P = ±1) if all vortices have 
positive/negative sign. FIG. shows the time evolution 
of the total number of vortices A^tot(^) = and 

of the polarization P{t) under influence of the obstacles 
(present throughout the whole evolution) with angular 
velocity uj. The top part (a) of the figure shows that 
A^tot(^) increases with uj. It is apparent that, by choosing 
cj, we can control the polarization. We use this tunable 
mechanism to create initial vortex distribution (without 
the pins) as shown in FIG. which plots Atot(^) and 
P{t) for initial states taken from instant t = 82 of FIG.|^ 
Clearly, by tuning uj we can produce a condensate free 


of external holes (the giant vortex or the obstacles) with 
approximately the desired vortex polarization. 



FIG. 3. (Golor online), (a) total number of vortices Wot; 
(b) polarization P vs time t, from initial states created at t = 
82 in the previous stirring process (FIG.[^, labeled by the an¬ 
gular velocities which generated them. The pins are removed 
and we evolve those states longer in time to study the vortex 
number decay. In (a), at t ^ 150, from bottom to top, the 
curves refer respectively to Ct; = 0, 7r/4, tt/IG, tt/G, and tt/S. In 
(b), from top to bottom, the curves are labeled as in FIG. 2 
(b). 

By numerically detecting each vortex and its trajec¬ 
tory, we determine Wot at each time step. By subtraction 
from the initial total vortex number, Nq = Wtot(O), we 
can infer the number of vortices which have drifted out of 
the condensate, W^r, and the number of vortices which 
have disappeared in annihilation events, Wan, colliding 
with vortices of opposite sign. We find that such vortex- 
antivortex annihilation events generate density waves, as 
already reported [Hiiaiia], turning kinetic energy into 
sound energy. The reverse mechanism is also possible m 
and in our 2D case takes the form of vortex-antivortex 
creation events, which we observe. Creation events oc¬ 
cur when the motion of the vortices induces a sufficiently 
deep density wave, or when a large amplitude wave ap¬ 
proaches the edge of the condensate where the local speed 
of sound c is less than in the central region. We have 
also observed annihilations events immediately followed 
by creation events: this sequence happens when a vor¬ 
tex collides with an antivortex, producing a large sound 
wave, which almost immediately generates a new vortex- 
antivortex pair, due to the changing value of the local ra¬ 
tio Vc/c; this effect happens near the condensate’s edge. 
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FIG. 4. (Color online) Positive and negative vortex number N~^ and N~ decay as a function of time t for cases (a) a; = 0; (b) 
u = 7r/16; (c) uj = tt/S; (d) u = tt/G; (e) u = 7r/4; and (f) phase-imprinted P = 0. The dashed lines are the respective fits for 
the numerical data (full lines) with the htting parameters Fi and F 2 appearing at top part of each plot. 


B. Polarized Turbulence Decay 

Starting from t > 82 (when we remove the pins and 
start a new simulation) we can examine whether there 
is a simple law for turbulence decay in 2D condensates. 
It has been suggested usmi] that the decay rate of the 
total number of vortices is not exponential and can be 
phenomenologically described by the logistic equation 

^ = -riVtot - (2) 

where the linear term should refer to vortex drifting out 
of the condensate, the non-linear term arises from vortex- 
antivortex annihilation events, and the coefficients Fi 
and F 2 are rates to be determined. We find that the 
solution of the logistic equation fits our decays for t > S2 
(after pins removal) fairly well. However, the fitting pa¬ 
rameters for the linear rate F 1 were shown to be negative 
in most cases, corresponding to positive growth. After 
the pins are removed, there is no reason or evidence to 
expect a vorticity source term, apart from occasional cre¬ 
ation of vortex-antivortex pairs mentioned above. There¬ 
fore, a naive association of the linear term of the logistic 
equation with vortex drifting out of the condensate due 
to dissipation effects is not appropriate in our case. 

Instead we propose a modified approach to the phe¬ 
nomenological description, focusing on changing the non¬ 
linear contribution to the rate equation. It is clear that 
we must incorporate polarization in the description. In 


Ref. [16], the authors include the non-linear term based 
on the intuitive assumption that annihilation processes 
are likely to occur with a rough dependence on the num¬ 
ber of dipole pairs (composed of a positive and a negative 
vortex) which can be formed. In a zero-polarized system 
this is of order oc N‘^. Following a similar reasoning, for 
a polarized system, the non-linear term should then be 
oc N~^N~. This consideration takes into account the fact 
that fewer vortex-dipole pairs are formed if P 7 ^ 0 and 
implicitly allows for polarization’s time dependence. 

With this motivation, we found that the phenomeno¬ 
logical model based on the coupled differential equations 


= -ri(v+)3/2 _ r2iN+N-f, 

_ (3) 

= -ri(V-)3/2 - T2{N+N-f, 

fitted well our simulation curves, giving non-negative val¬ 
ues for both fitting parameters Fi and F 2 , as show the 
respective values inserted in the plots of FIG. 

The 3/2 scaling for the first terms on the right-hand 
side of Eqs. arises from a simple argument which as¬ 
sumes that: (1) a singly-charged vortex (with either neg¬ 
ative or positive circulation, k = ^1) moves radially with 
a typical velocity proportional to Vr = l/(27r/), where 
I = l/\/n is the inter-vortex spacing and n = 
is the vortex density in a condensate of radius P; (2) 
vortices are expelled from the cloud in a proportion P, 
which in first approximation is independent of the vortex 
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number N but may depend on the dissipative parame¬ 
ter 7 . Therefore, in an infinitesimally thin annulus of 
area dA = 2'KRvj.dt close to the border of the condensate 
there is a loss of dN = —SndA, and we get the scaling 
dN/dt (X A^3/2. 

In order to account for the initial steep descent of some 
of the curves, it was necessary to square the polarization- 
dependent term of the equation, i.e. otherwise 

the fits have shown to be very poor. Therefore, summing- 
up Eqs. the total number of vortices decays non- 
trivially as a function of polarization P according to 

^ = -ri77 - r2/(i)iv7, (4) 

where the time-dependent polarization appears in the 
function f{t) = (P(t)^ —1)^/8. The same scaling for 
the vortex number decay rate was justified heuristically 
in the context of a quenched 2D homogeneous system in 
Ref. [45]. 

From FIG. we can identify cases (b) and (e) (cj = 
7r/16 and u; = 7r/4), (c) and (d) {uj = ir/S and uj = tt/G), 
respectively as opposite-polarization counterparts; their 
decay curves are similar in behavior, while having a 
polarization mirror-symmetry. The main difference be¬ 
tween these curves and their counterparts is the steeper 
number decay in initial times for (c) and (e). Calcu¬ 
lating the number of vortices decayed due to annihila¬ 
tions, we found that this loss was always considerably 
less than the drifting. The latter is strongly induced by 
vortices interactions, which push each other out of the 
cloud. Since our chosen dissipation was shown to be inef¬ 
fective in making vortices spiral-out of the condensate in 
the time-scale studied, we conclude that both linear and 
non-linear terms in Eq. 0 have origins on vortex inter¬ 
actions. Case (e) {uj = tt/A) illustrates well this need for 
steeper than quadratic term in the rate equation, since 
it characterizes a purely non-linear decay (i.e. Ti = 0 ). 

Finally, in order to compare models from Eq. ([^ with 
Eq. 0 we performed numerical experiments in which, 
rather than creating vorticity with the giant vortex-pins 
set up here proposed, we simply numerically imprint a 
given initial number of vortices uniformly at random po¬ 
sition onto the same harmonically trapped condensate. 
We obtain essentially the following results: the polariza¬ 
tion approximately retains its initial value P = 0 (see 
plot ^ in FIG. 1^, and a reasonable fit is obtained using 
Eq. Q. As opposed to the polarized cases, we find non¬ 
negative rates. However, we see in the comparison shown 
in FIG.j^that Eq. Q fits the curve better, showing that 
in this case the non-linear term should be oc instead 
of oc N‘^ to account for the steeper decay. Similarly to the 
polarized cases, drifting has shown to be the main mech¬ 
anism of vortex loss. Eq. 0 has fitted well the cases 
studied in uniiiH], which differs to our P = 0 case not 
in polarization but rather in the initial total number of 
vortices (in their case ^ 60 as opposed to ^ 140 in ours). 
Therefore, we attribute the strong non-linear effects and 
the departure from the quadratic nature (which is con- 



FIG. 5. (Color online) Total number of vortices, Ntot vs time 
t for the case where an unpolarized P = 0 vortex distribution 
was created through phase-imprinting, in order to compare 
fits given by Eq. § (black dot-dashed line) and Eq. 0 (blue- 
dashed line) 


sistent with the ‘ultra-quantum’ decay observed HE] in 
superfiuid helium, where system’s finiteness is not an is¬ 
sue) to vortex mutual interaction in a system which is 
finite. The finiteness imposes a certain limit to the num¬ 
ber of vortices which can be accommodated in the cloud. 

In summary, for both polarized case and the particu¬ 
lar unpolarized case where the vortex density is high, our 
proposed modified decay rate equation, Eq. 0 success¬ 
fully describes the total vortex number evolution. Fur¬ 
thermore, it adequately applies to our particular dissipa¬ 
tive simulations. 


V. CONCLUSION 

We have presented a new scheme for generating 2 D 
quantum turbulence in atomic condensates which allows 
control over the polarization of the flow, equivalent to 
the net rotation of a turbulent ordinary fluid. Using this 
experimentally feasible scheme, we have examined the de¬ 
cay of the turbulence and the vortex interactions (vortex- 
antivortex creation and annihilation) which take place 
in the condensate. These results modeling the decay of 
the number of vortices using a new model equation that 
takes into account polarization and its time-dependence; 
we have found that, in this context, it offers more rea¬ 
sonable fits than the logistic equation proposed by m- 
The reason being it includes polarization in the descrip¬ 
tion and also accounts for more non-linear (steeper than 
quadratic) effects that have proven to be important in 
the initial decay process. 

In short, the inclusion of the time-dependent polariza- 
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tion parameter to the rate equation of the vortex number 
seems to capture the complete dynamics. It can offer ex¬ 
perimentalists a hint on the ratio of positive and negative 
vortices in a 2D system, which is still an experimental 
challenge nowadays. 
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